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We present a simple, relatively efficient method to predict the nonvolatility 
of dual-dielectric, charge-storage cells. Using this method, charge-retention 
times of several hundred years can be predicted unambiguously from ex- 
periments of several days' duration made at elevated temperatures and with 
externally applied, accelerating fields. The method is first presented in the 
context of a simple, physically reasonable model of the bias and temperature 
dependence of the characteristic relaxation time of the device. It is then 
used to analyse a particular device structure. At 80°C, we predict that this 
device can maintain a flatband-voltage shift in excess of 4 volts for approxi- 
mately 500 years. Our analysis suggests that this method can be applied 
to a variety of dual-dielectric, charge-storage cells to predict their non- 
volatilities. 

I. INTRODUCTION 

The retention time 1-9 of charge in dual-dielectric, charge-storage 
cells 10-12 (DDC's) is of central importance in evaluating and comparing 
these devices. An ideal device would hold its charge at a constant level 
indefinitely. While, of course, such an ideal device is physically im- 
possible, nonideal, charge-storage memory cells have been fabricated 
that are expected to have minimum charge-retention times on the 
order of tens of years. And, as such devices are improved, even longer 
retention times can be expected. The question that naturally arises is 
how to evaluate such devices in a few days to predict their charge- 
retention times. 

The primary purpose of this paper is to describe a method of bias- 
temperature stressing in which the decay of the stored charge is greatly 
enhanced and, as a result, from which one can predict the charge- 
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retention time of the device under normal operating conditions (zero 
bias and room, or somewhat higher, temperature). We stress at the 
outset, moreover, that our primary concern is charge retentivity and 
not device reliability. Bias-temperature-stress methods are often used 
to good advantage in reliability studies to artificially reduce the 
lifetime of the device. Here, we are using similar methods to arti- 
ficially reduce the retention time of charge stored on the device. 
For reliability studies, one makes use of models in which the device 
lifetime is governed by temperature and applied bias. Owing to the 
complexity of the aging, however, these models are often quite tenta- 
tive. Here we use a well-defined model in which the rate of decay of the 
stored charge is governed by temperature and bias. As in reliability 
studies, experimental results are used to determine the parameters of 
the model. While the bias dependence of the rate of decay of stored 
charge has been noted previously, 7 the work we report on here is the 
first in which bias and temperature stressing have been used simultane- 
ously to predict charge retention under normal operating conditions. 
We feel obligated, therefore, to elaborate our approach and findings in 
some detail. 

At the heart of our method of prediction are the observations (i) 
that there exists an "envelope" function that sets an upper bound on 
the total stored charge that can be present in the device at any given 
time, independent of initial conditions, (it) that this envelope function 
is determined primarily by the characteristic relaxation of the device, 
and (m) that this relaxation is a strong function of temperature and 
electric field. By focusing attention on the envelope function, which 
properly indicates the long-time decay of the stored charge, we avoid 
incorrect predictions of decay time based on extrapolations of the 
initial portion of the charge-decay curve. (If the charge decay is 
plotted versus log time, such predictions are overoptimistic ; if it is 
plotted versus time, such predictions are overpessimistic.) From mea- 
sured values of the characteristic relaxation obtained at elevated tem- 
peratures and under applied bias, relaxation times of interest at room 
temperature and under zero bias can be predicted. Owing to the non- 
linearities inherent in the charge decay, these zero-bias, room-tem- 
perature relaxation times are functions of the initial stored charge; 
the larger the initial stored charge, the smaller the relaxation time. 
For example, for one version of our device 13 we find that, for an initial 
charge corresponding to a flatband voltage of 10 V, a relaxation time 
of 3 -10* years is predicted; for 7 V, we predict 6 10 4 years. Using our 
results, predictions of relaxation times can also be made for devices 
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operated under bias and at other than room temperature. For ex- 
ample, at 80 °C under zero bias we predict for 10 V a relaxation time 
of 100 years and for 7 V 300 years. Linear extrapolations of the initial 
portion (t « r n ux) of the charge-decay curve would result in misleading 
estimates of charge-retention times on the order of 10 20 years. 

It is important at the outset to stress that we are concerned in this 
paper with devices 10 with "thick" oxide layers; that is, with oxide 
layers between the semiconductor and the charge-storage sites suffi- 
ciently thick that the rate of decay of charge via back-tunneling 
through the oxide layer is small compared to the rate of decay of 
charge through the insulator layer between the storage sites and the 
gate. While one can accelerate the back-tunneling charge decay of thin- 
oxide devices 7 by applying an electric field, we believe that, unless one 
has an exceptionally good means of characterizing the specific tunneling 
processes of interest, such information is of little direct value in esti- 
mating the decay in the absence of such field, i.e., under normal oper- 
ating conditions. The applied voltage, of course, can accelerate or 
decelerate decay through the oxide layer or the insulator layer, depend- 
ing upon its polarity. If electrons are stored, a positive gate voltage 
enhances the decay in thick-oxide devices because the decay is through 
the insulator, whereas the same voltage reduces the decay in thin-oxide 
devices 7 where the primary decay is back through the oxide. One 
reason for studying thick-oxide devices is to determine the limits on 
the retention time of charge-storage devices imposed by charge decay 
through the insulator layer. Such limits are of considerable interest, 
especially for devices whose back-tunneling decay is sufficiently low 
that, based on this decay mechanism alone, one might exaggerate the 
device's charge-storage capabilities. 

In this paper we first discuss in general terms the physical processes 
that lead to the decay of the stored charge. A relatively thick oxide 
layer is used so that the primary discharging current is through the 
insulator layer. This current through the insulator is found to be char- 
acterizable as a thermally activated flow of charge via defects with a 
very low, but strongly field-dependent, mobility. We then discuss how 
this strong temperature and electric-field dependence of the decay 
current can be used to controllably accelerate the discharging of the 
DDC. A simple, physically reasonable, empirical model is introduced 
to explain the experimentally measured decay of the flatband-voltage 
shift in time with temperature and externally applied bias as param- 
eters. The mathematical expression for the decay current is sufficiently 
simple that all quantities of interest, particularly device relaxation 
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times, can be calculated analytically. From the high-temperature, 
high-accelerating-field results, it is possible to predict low-temperature, 
zero-bias behavior. We can thus unambiguously determine nonvola- 
tility on the order of tens or hundreds of years on the basis of experi- 
ments performed in several days or weeks. 

Before proceeding, it must be emphasized that the method of bias- 
temperature stressing that we develop here is not limited to devices 
which behave according to our simple empirical model. Indeed, as we 
point out in Appendix A, the central ideas at the heart of the method 
are valid, within certain limits, to a variety of models : the method is, 
within certain bounds, insensitive to the details of the actual decay 
process. We discuss the physical processes that we believe are opera- 
tive in our devices first, however, because this will permit the reader 
to familiarize himself with the actual devices to which our method is 
then applied. Although our results indicate the general features of a 
particular, detailed decay model, we emphasize that we are not pri- 
marily concerned with establishing the existence of such a model. 

II. PHYSICAL PROCESSES 

The structure and fabrication of the memory devices studied here 
have been described in detail elsewhere. 13 For our purposes here, it 
suffices to note that a typical device consists of the following layers 
(see Fig. 1) : a metallic contact (taken fixed at zero voltage), n- (or p-) 
type silicon, 70 to 200 A of Si0 2 (referred to as the oxide layer, or 
simply the oxide), a set of suitable dopant-induced storage states 
(referred to as the storage states), 400 to 500 A of Al 2 O s or Si 3 N 4 
(referred to as the insulating layer, or simply the insulator), and a 
metallic contact (referred to as the gate). Inclusion of the specific 
states, which are the interfacial dopant-induced states, makes it 
possible to store either electrons or holes, whichever is desired, as well 
as to provide very well-defined storage sites for the elementary charges. 
The oxide and the insulator are sufficiently thick that tunneling from 
one side of either to the other side can be neglected. 

The device can be charged either negative or positive as follows. 
Negative charge can be stored by driving the gate to a high positive 
voltage (?^50 V) for a short period of time (« 100 ms). Some electrons 
in the Si tunnel into the conduction band of the oxide, traverse the 
oxide, and then are trapped in the storage states or pass completely 
through the device to the gate electrode. Positive charge is stored by 
driving the gate negative to force electrons out of the storage states, 
through the oxide, and into the semiconductor. 
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Fig. 1 — Schematic of the structure of the MIOS memory device under five differ- 
ent, applied-bias conditions. The location of the interfacial, dopant-induced, charge- 
storage states is indicated. They are assumed to be charged negative. 

Our experimental results indicate that these multilayered memory 
devices discharge as follows. (Refer to Fig. 1 in which a device is shown 
under five different biases. Note, in particular, Fig. la, in which the 
semiconductor is inverted near the Si0 2 .) The (negative) stored charge 
produces large electric fields in both the oxide and the insulator. In the 
first few hundred milliseconds, some fraction of the charge is lost 
through thermionic emission or tunneling into the conduction bands of 
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the oxide and/or the insulator. For most of the stored charge, however, 
the trapping levels are sufficiently deep that these processes are very 
improbable. Charge can tunnel into traps in the insulator, and through 
such imperfections a discharging current can pass into the conduction 
band of the insulator and flow to the gate. Some charge can also pass 
completely through the insulator in the forbidden band by means of 
these traps. In addition, in a similar manner, a discharging current can 
be present in the oxide and flow to the Si. 

The size of the discharging current is expected to be strongly field 
dependent. 1-9 In addition, we note that increasing the (positive) gate 
voltage enhances the electric field in the insulator and decreases this 
field in the oxide. The manyfold increase in the decay current which 
occurs when a positive voltage is applied to the gate thus implies that 
the primary discharging current is passing through the insulator. In- 
versely, when a negative voltage is applied to the gate, a decrease in 
the decay current is observed. In this case, the change in the applied 
fields is such as to increase the current through the oxide and decrease 
the current through the insulator. A detailed study of the reverse gate- 
bias decay was not undertaken because of the long time intervals 
required even at high temperatures. For the present, it is sufficient to 
note that the total decay current increases with positive voltage and 
decreases with negative voltage applied to the gate. 

The physical details of the transport of the charge within the insu- 
lator can be narrowed down as follows. The strong temperature en- 
hancement of the decay of the charge indicates that thermal activation 
rather than tunneling plays the dominant role in enabling charge to be 
transferred toward the gate. The question of whether the charge passes 
through the insulator to the gate entirely within the forbidden band, 
or partially in the conduction band, is more difficult. The relatively 
low thermal activation energies observed (« 0.6 eV) suggest that the 
rate limiting portion of the transport is not associated with the con- 
duction band. If it were, much higher decay currents would be ex- 
pected. We infer, therefore, that the current is controlled by the traps 
in the bulk of the insulator. 

For charge transport from one trap to the next, one might expect 
the current to be proportional to exp [(E„ + qSa/2)/kT], according 
to Poole's law. 14 (Here E is an ionization energy, 6 is the electric field 
in the insulator, and a is the intertrap spacing.) However, Poole's 
law, or its extension and modification by Frenkel 15 ' 16 and others 17-19 
include only the thermal excitation of carriers and do not include 
velocity-versus-field effects connected with the transport from one site 
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to the next. 20 Current density is proportional to qnv, where n is the 
number of charges per unit volume and v is the carrier velocity. One 
expects n to be thermally activated; however, in an insulator, v can 
be a rapidly varying function of applied fields, 21-23 whether electron 
transport is in the conduction band or in a defect "band," as seems to 
be the case here. Thus, whereas for pure Frenkel-Poole behavior the 
excitation of electrons is the rate-limiting, field- and temperature- 
modulated process, for our devices both n and v are so modulated. 

Our experimental results indicate that the discharging current is 
proportional to the empirical expression 

8 {T) exp {-EJkT) exp [+«/«.(Dl (1) 

where S (T) is a very slowly varying, decreasing function of the tem- 
perature T. To the accuracy of the experiments, the activation energy 
E a is independent of the applied field. The decrease of S„(T) with 
temperature could be due to an increased overlap between polarized 
electronic states because of increased lattice vibrations. The electric 
field dependence can in general be expected to be more complicated 
than that given in eq. (1). Fortunately, eq. (1) suffices for our purposes, 
partially because the decay curves are relatively independent of the 
detailed dependence of the discharging current on electric field. See 
Appendix A for a discussion of more complicated field dependences. 

One final point should be made regarding the decay current. As part 
of the assumption that the rate-limiting process controlling the charge 
decay is the trap-controlled current within the bulk of the insulator, 
we have assumed that the number of carriers in transit is independent 
of time. In other words, as a trap in close proximity to the storage 
states loses its charge to a trap somewhat closer to the gate, the emptied 
trap is quickly refilled with a charge from a storage state. Thus, the 
average time for a charge to pass from a storage state to an empty trap 
is much less than the average time for a charge to pass from one trap 
to the next. This results from the high density of storage states, and 
makes the effective number of carriers available for transit independent 
of the magnitude of the stored charge and hence independent of time. 
This is reasonable throughout most of the decay owing to the large 
number of stored electrons. Near the end of the decay, however, when 
the amount of stored charge is much less than its initial value, the rate- 
limiting process may become the transfer of charge from the storage 
states into the insulator. Wc ignore this effect, since, by the time it 
becomes important, the fundamental relaxation time of the device will 
be well determined. 
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From the foregoing discussion, it is evident that we can artificially 
enhance the decay greatly by increasing the temperature and the 
electric field in the insulator. In so doing, however, we must be very 
careful that we do not excessively enhance charge-transport mecha- 
nisms that under normal operating conditions would play no sig- 
nificant role. For example, by applying too large a potential at the 
gate, the field in the oxide will be enhanced to such an extent that the 
storage layer will be slowly charged by electrons tunneling from the Si 
into the conduction band of the oxide. Excessive temperatures, on the 
other hand, may enhance decay into the insulator's conduction band, 
exaggerating this mode of charge transport. These effects are often 
easily detected experimentally, and they give rise to upper limits on 
the increases in temperature and electric field possible to enhance 
charge decay for measurement purposes. 

Although the decay of charge from the storage states can be en- 
hanced by increasing either the temperature or the insulator electric 
field, the enhancement from either is insufficient to bring the decay 
times into the range of days. However, this can be achieved by com- 
bining higher temperature with higher fields. Of course, we again must 
be careful to avoid introducing decay processes that under normal con- 
ditions would be insignificant. (An example of this would be Schottky 
emission.) Nonetheless, we find that we can obtain the entire charge 
decay curve within a few days (at most) without introducing extrane- 
ous decay mechanisms. From these curves, we can predict the normal 
decay curve and hence the charge-retention time of the memory ele- 
ment. In this way we avoid having to determine this retention time by 
extrapolating the initial portion of the time dependence of the charge 
decay. In the next section we see how the nonlinear dependence of the 
rate of decay on the quantity of stored charge leads to characteristic, 
charge-decay curves from which one can predict the nonvolatility of 
DDC's. 

In those cases where certain alternative conduction mechanisms arc 
activated, this method may or may not work. In some devices, the 
charge decay is via Fowler-Nordheim tunneling from the storage states 
into the conduction band of the oxide or the insulator, or it is via direct 
tunneling from these states into the semiconductor. Neither of these 
decay currents will be affected by changes in the temperature; how- 
ever, both such currents will be strongly modulated by applied bias. 
In the case of Fowler-Nordheim tunneling, our method is directly 
applicable. For direct tunneling, however, in which the primary decay 
is via Si-Si0 2 interfacial states, we must be very careful as we increase 
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the bias not to enhance the tunneling into semiconductor states in the 
valence or conduction band, which under normal conditions would 
play no important role in the discharging of the DDC. Perhaps one can 
modulate the decay by altering the thickness through which the elec- 
trons tunnel, either by physically squeezing the sample or by con- 
structing samples of varying thicknesses. 

III. EMPIRICAL MODEL AND MEASUREMENT TECHNIQUE 

Let us assume that we have some (negative) charge stored in the 
interfacial, dopant-induced states. The quantity of charge stored in 
these states is proportional to the measurable flatband voltage V fb . Also 
proportional to Vfb is the field in the insulator under zero-bias con- 
ditions. In the presence of a finite bias voltage Vb, an additional field 
proportional to the bias voltage will be linearly superimposed on the 
zero-bias field (with a different constant of proportionality). Thus, 
using expression (1) and expressing fields in terms of voltages, we 
write for the time rate of change of V fb) the flatband voltage, 



dVfb = V.jT) 
dt t(T) 



f V fb + bV b 1 
L V.(T) J ' 



**P ^H^ > (2) 



where V, and r are experimentally determined functions of temperature 
only and b is a relative constant of proportionality (sec Appendix B). 
Equation (2) presumes decay via the insulator ; for decay through the 
oxide, one would have the factor (V fb — V b ) in place of the factor 
(V/b + bVb). This empirically satisfactory expression is reasonable 
physically, as explained in Section II. It also permits a simple, analytic 
treatment of all significant features of the decay. Should the flatband 
voltage decay according to a relation other than (2), one can still 
employ bias-temperature stressing. This is indicated in Appendix A. 
Before solving eq. (2), we should clarify two points, (i) We do not 
require eq. (2) to be valid for t < 1 minute. It may be valid there but, 
being interested in the long-time decay of the charge, all we need to do 
is to integrate eq. (2) from some time h, on the order of minutes after 
charging, when the experiments are begun, (ii) Since dV /b /dt does not 
vanish as (V /b + bV b ) ->0, eq. (2) is clearly not valid as t — > «>. As 
remarked in Section II, this difficulty results from assuming that there 
are always a large number of stored charges, whereas in fact near the 
end of the decay this number becomes small. Again, this difficulty need 
not concern us, since we do not have to consider the leakage after the 
bulk of the stored charge has decayed away. (See Appendix A for a 
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more general treatment.) With these limitations in mind, we proceed 
to solve eq. (2). 

The solution to eq. (2) follows at once after a separation of variables 
and an integration from (t h Vi) to (t, V fb ). The resulting V fb (t) for the 
decay of the flatband voltage as a function of time t at temperature T 
and applied bias V b is given by 

V, b (t) = 7,(r)log e {(exp [-Vi/V.Cr)] + [« - ti)/r(T)2 

■explbV b /V 8 (T)2)- 1 }. (3) 

For a single value of T, this function is sketched in Fig. 2 as a function 
of log™ (t/ti) for several V b and for hypothetical, illustrative values of 
V„ t, and b. The purpose of plotting as a function of log™ (t/h) is to 
call attention to the actual long-time behavior of V fb (t). For h < t 
« t(T) exp [- (Vi + bVt)/V.(T)2, V fb (t) is relatively flat. However, 
for t » t(T) exp [- (Vi + bV b )/V.(T)] t V fb (t) is given approximately 
by V,(t) where, if h is any convenient scale factor, 

V.(t) = V.(T) log, [r(T)A'i] - bV b - V,{T) log e it/h). (4) 

The function V e (t) is the previously mentioned envelope function for 
each T and V b . It is shown dashed in Fig. 2. (Note log e x = logio x/ 
logio e.) Although h may be assigned the same value as t\, the two 




102 
TIME t IN MINUTES 



Fig. 2 — Schematic illustrating a possible decay of the flatband voltage V/h(jt) as a 
function of logio (t/k) in a hypothetical device for four different applied biases V b 
at a fixed temperature T. For simplicity, we set b = 1. Also shown are the envelope 
functions V c (t) associated with each Vb- 
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must be distinguished. Our h refers to an initial condition, while V,(t) 
is independent of the initial conditions, t\ and Vi, and satisfies V f b(t) 
< 7,(0 (see Appendix C) and 

V fb (t) « 7.(0, t » r exp [- (7, 4- bV b )/V a J (5) 

When 7.(0 is plotted versus log™ (t/h), V, is the slope and r is the 
intercept [F«(r) = 0] for zero bias (V b = 0). From the experimental 
V f b (t) curves obtained for several Vb at a given T, it is possible to 
readily determine t(T), V,(T), and b, the three parameters in eq. (4). 
From the envelope function V e {t) drawn tangent to V/b(t) as shown in 
Fig. 2, we obtain V,(T) from the slope of V e (t), b from the difference 
between V„(t) (for fixed t) at different Vb, and t(T) from the value t 2 
of t where 7.(0 = 0: 

log e [r(T)/U = log. ( |) + ^y (6) 

In Fig. 2, t 2 for 7& = 8 V is shown (fe = 5-10+ 3 minutes). It is a test 
of the form of eq. (2) that V B (T), b, and t(T) be independent of V b 
to within experimental accuracy. That this is indeed the case is dis- 
cussed in Section IV. We now repeat the above for other temperatures 
in the range 150°C ^ T ^ 300°C. We find that both r(T) and V,(T) 
appear to be "thermally activated" : 

r(T) = t exp (Ea/kT), V.(T) = V ao exp (E./kT). (7) 

Since the slope V e may arise from a combination of physical processes, 
identifying it with a thermally activated process may be somewhat 
misleading. (The parameter b is a weakly decreasing function of T. 
Since we are interested only in predicting zero bias (Vb = 0) behavior, 
we need not concern ourselves further with extrapolating b to room 
temperature. In other words, we can for each temperature use the b 
measured at that temperature to extrapolate to zero bias. Once we 
know the zero-bias result, we can extrapolate to the desired tempera- 
ture. Since b enters only as bVb and Vb is zero, this latter extrapolation 
can be performed without knowledge of b.) The quantities r , E a , V t0 , 
and E, are obtained from the usual (1/T) plots. We assume, and 
indeed from our discussion in Section II it is not unrealistic to expect, 
that we may extrapolate t and V, to room temperature using (7). We 
must now relate these quantities to the characteristic relaxation time 
t relax, the charge retention time, of the storage device. 

In Fig. 3 we define r re iax(T r ) graphically, again for hypothetical 
device parameters. It is simply the "roll-off" time of 7/j(0 : 7,(T re i a x) 
= 7i, where V\ = V /b (ti). [As shown in Appendix C, Vi « V fb (0) 
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Fig. 3 — Schematic illustrating a possible decay of the flatband voltage V(t) in a 
hypothetical device for room temperature (T = T ) and zero applied bias (Vb = 0) 
for two different quantities of stored charge, Fi and V[. The size of the characteristic 
relaxation time T re i«i is indicated for both cases. Note that V x — V\ = V l — V° 
= 0.69 V,. 

for h « Treux.] Solving (4) for T re iax at room temperature T for V b => 0, 
we obtain 

Tre.ax = t(T ) eX P [- V t /V. (T D )J (8) 

Indicative of the nonlinearities inherent in the charge-decay process, 
Trelax is a function of the initial stored charge, and hence the initial 
flatband voltage V h increasing as V\ decreases. (Compare T rc i B x and 
rreiax corresponding to different amounts of initial stored charge Vi 
and Vi, respectively, in Fig. 3.) Again we emphasize that T re i ax cannot 
be obtained from (linearly) extrapolating V(t) based on its behavior 
for t « T re iax- Since T re iax is on the order of 10 4 years at room tempera- 
ture for devices reported on here, it is doubtful whether room-tem- 
perature, zero-bias experiments alone will be able to predict Treiax- 

Although Treiax provides a measure of the charge retention time of 
the device, we must, in addition, know how far V S h has decayed below 
Vi by the time T re i ax before full significance can be attached to this 
definition of characteristic decay time. Returning to eq. (3) and letting 
t x = 0, Vi = V', V b = 0, we obtain for V fb {t) under zero bias 



V 



n 



(0 = F .log.[^-(l + ^;)- 1 ] 



(9) 
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In arriving at eq. (9), we have used eq. (8). Thus, 

7/*(r„u«) = V\- 7, log, 2 

= VI- 0.697,. (10) 

We thus obtain the very important result that, at the roll-off time 
Treiax, V f0 has dropped 0.697, volts below its initial value. This amount 
is independent of V°, the initial flatband voltage. Of course, 7 /6 (r re iax) 
must be several volts positive so that eq. (2) will still be valid. In addi- 
tion, 0.697, thus defines the decay that must be taken into account 
when including this device in circuits. In determining 7, from the 
envelope of a decay curve, one must be careful to convert from log™ 
to log,: 

T. = ™-™> X 0.434, 
logio (h/h) 

where 0.434 = logio e. 

Another characteristic decay time is the time at which the flatband 
voltage drops below a certain margin voltage V m . We refer to this time 
as the margin time, T margin , or simply r m . It follows at once from eq. (9) 
that r m as a function of 7 m and V\ is given by 

r m = r(T){exp [-7 m /7,(D] - exp [- 7 1 /7,CT)]| (11) 

as a function of T under zero bias conditions. A plot of T ma r K in is pre- 
sented in the next section. 

IV. EXPERIMENTAL RESULTS 

Bias-temperature-enhanced, charge-relaxation experiments were 
carried out on the recently devised dual-dielectric, charge-storage 
cells. 13 Aluminum oxide (Al 2 Oa, e = 9) was used as the insulator. The 
thickness d of the Si0 2 (« = 4) was 125 A, and that of the A1 2 3 , d i} 
550 A. Thus, b = C /d [see eq. (25)] is predicted to be given by 

b _ Co _ esio^l d, = 20 
Ci d e A i,o 3 ^ 

Experimental values for b ranging from 0.8 to 1.2 were obtained over 
the temperature range of 150°C ^T ^ 300 °C. The reason for this 
discrepancy is not known. Some possible explanations are discussed in 
Appendix D. The area of the devices was about 1 • 10 5 /xm 2 , and the 
doping of the n-type silicon was about 6 • 10 15 per cm 3 . Tungsten of a 
density of about 2 to 3 • 10 14 per cm 2 was used to produce the interfacial 
storage sites. 
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Table I 



Sample* 


T 

(°0) 


i/r-io» 


v b 

(volts) 


v. 

(volts) 


(minutes) 


b 


r(T) 
(minutes) 


H15-4 
H15-1 
H15-1 
H15-4 
H15-4 
H15-1 
H15-1 


150 
200 

250 

300 


2.36 
2.11 

1.91 

1.74 


15 

{It 

u 


2.44 
1.93 
1.92 
1.76 
1.76 
1.53 
1.52 


2.6-10* 

i.3 ion 

7.7-10*/ 
2.7-10*1 
4.7 Wj 

4.6 -ion 

3.3 10V 


1.2+ 
1.12 

1.03 

0.8 


3.75 10' 
4.3 -10 8 

1.7 10« 

6.2 10 6 



* H15 is the wafer label and 1 or 4 is the chip label. 



Table I presents a summary of the experimental results which we 
have analyzed most carefully. Four representative charge decay curves 
are plotted in Fig. 4. The envelope function V e (t) associated with 
each curve is shown as well. From the V e (t) we obtain the data given 
in Table I. One might do better by attempting to fit the actual experi- 
mental curves with eq. (3), but we have found V e (t) adequate for our 
purposes. 

In studying Fig. 4 one may wonder why the slope of the decay curves 
seem to diminish at higher temperatures. In fact, as is clear from eqs. 




102 103 

TIME t IN MINUTES 



Fig. 4 — Experimental results under four different conditions: T = 150°C, 
V b = 15 V; 250°C, 7 V; 250°C, 10 V; 300°C, 10 V. In each case, the element was 
charged by setting Vb = 50 V for 100 its. 
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(2) and (7), \dV fb /dt\ increases as T increases, as one would expect. 
One must remember that we are plotting V fb (t) versus log™ (t/ti). 
Thus, the enhanced decay has shifted V fb (t) to lower t, reducing the 
apparent decay of V fb when plotted on a log™ t scale. The authors 
must admit to having expended some effort themselves getting used to 
plotting V Jb versus log t rather than log V /b versus t, as is common in 
many decay problems. 

By plotting the values of t and V s given in Table I as functions of 
1/T, the activation energies E a = 0.61 eV for r and E, = 0.065 eV for 
V, can be determined [see Fig. 5 and eq. (7)]. From Fig. 5 we can 
also calculate r and V so [eq. (7)] : r = 2.0 minutes and V, = 0.41 V. 
Using the E a , t , E„ V, so determined, we use eq. (7) to extrapolate 
t(T ) and V,(T„), where T is room temperature (290°K, 17 °C). We 
find t(T ) = 7.7 X 10 10 minutes and V.(T ) = 5.5 V. With these 
values, we can predict the room-temperature, zero-bias relaxation 
times using eq. (8). 

The activation energy E a of 0.61 eV is comparable to activation 
energies ranging from 0.4 to 1.2 eV reported for Frenkel-Poole con- 
duction through oxide layers. 24 " 28 We conclude, therefore, that E a may 
be interpreted as being associated with an activation process. Although 
E. is relatively small, being on the order of several kT's, its existence 
does result in an appreciable variation of V B with T (as observed experi- 
mentally). We do not have a simple independent quantitative explana- 
tion of its magnitude and caution the reader against interpreting E a as 
being associated with a simple, thermally activated process. 

As discussed in Section III, the charge-retention time of our dual- 
dielectric, charge-storage cell is well-characterized by the device 
"roll-off," or relaxation time, as defined in Fig. 3. Owing to the non- 
linear dependence inherent in the charge decay, r re i a x depends upon 
the "initially" stored charge, which we denote in terms of flatband 
voltages by 7, = V^ih). In Fig. 6 we plot r re i ax as a function of V x 
for several temperatures of interest. For room temperature and a Vi 
of 10 V, T r eux is 3 ■ 10 4 years ; for 7 V, 6 • 10 4 years. That r re i a x is a func- 
tion of stored charge should be carefully noted. Lastly, we note that, 
for room temperature at r»i«, V /b has dropped (0.69) (5.5) = 3.8 V 
below V t . At 80 °C the drop is only 2.4 V owing to the reduced value, 
of V. ( = 3.4 V at 80°C). However, now t„i« for 10 V is only 100 years, 
and for 7 V, 300 years. 

In Fig. 7 we plot r marg i n as a function of Vi for zero bias and for the 
same temperature used in Fig. 6. As is obvious physically, the larger 
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1.9 2.0 2.1 2.2 

RECIPROCAL TEMPERATURE (1/T-1Q3) 



Fig. 5 — Experimentally determined parameters t and V. plotted versus \/T to 
determine the activation energies E a and E.. 



the initial stored charge (initial flatband voltage Vi), the longer the 
time required for the flatband voltage to drop to the margin voltage 
V n , here taken to be 4 V. The horizontal lines give the upper limit or 
largest possible T marg i n for the given temperature. For V x between 4 and 
5 V, r m goes rapidly to zero, as is clear from Fig. 3. For V\ less than 
4 V, r m is obviously meaningless. 
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10H 



1013 - 



1012 - 



V. 10"- 



•• ioio _ 



109 



108 - 



10? 



w-Trne-vi"* 



T= 17°C 



(ROOM TEMPERATURE) 




4 6 8 10 

INITIAL FLATBAND VOLTAGE V, IN VOLTS 



108 



10' 



106 



-10* 



Fig. 6 — The predicted roll-off or device relaxation time for zero-applied bias 
(V b = 0) and temperatures of -40°C, 17°C, 80°C, and 100°C plotted as a function 
of "initial" flatband voltage Vi = V(U), where ti = 1 minute. 



V. CONCLUSION 

In this paper we have discussed a method of predicting the retention 
time of charge in dual-dielectric, charge-storage cells. The method is 
based on the realization that an envelope function V e (t) exists, that 
this function is determined primarily by the characteristic relaxation 
of the device, and that this relaxation is a strong function of tempera- 
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10" 



1013 



1012 



1/3 1Q11 



1010 
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108 
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W = T|Tl{a-WTU-V,/V s ITl} 
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(ROOM TEMPERATURE) 



/ 




4 6 8 

INITIAL FLATBAND VOLTAGE Vi IN VOLTS 



105 



104 



103 



102 



12 



Fig. 7 — The predicted margin time for zero-applied bias (Vb = 0) and tempera- 
tures of -40°C, 17°C, 80°C, and 100°C plotted as a function of "initial" flatband 
voltage Vi = V(0i where h = 1 minute. The horizontal line above each curve is 
the upper limit for T mor gi n as Vi is taken larger and larger. 

ture and applied bias. By choosing an appropriate function to represent 
the dependence of the discharging current on the applied bias and 
stored charge, it is possible to discuss all the interesting features of the 
time-dependent charge decay in terms of simple mathematical expres- 
sions. We were thereby able to predict, on the basis of experimental 
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data obtained under bias-temperature stressing of representative 
devices, the time dependence of the shift in the flatband voltage at 
zero bias and operational temperatures (room temperature to 80 °C). 
In particular, we can characterize the nonvolatility of the DDC's in 
terms of a roll-off or relaxation time r T and a margin time r m , both 
dependent on the level of initial charging. Although these times are on 
the order of 10 4 years, they were determined from experimental data 
taken over a period of only a few weeks. This was possible because our 
method of bias-temperature stressing greatly enhances the dominant 
discharging processes, speeding up the charge decay without intro- 
ducing significant decay from otherwise insignificant decay modes. 
Because of this, we feel that our method will greatly facilitate the 
realistic prediction of the nonvolatility of DDC's. 
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APPENDIX A 

Equations (1) and (2) in the text are admittedly oversimplified, 
despite the ease they render in interpreting the data. Nonetheless, the 
characteristic decay V{t) discussed in the text is approximately ob- 
tained in a variety of cases, e.g., tunneling or Frenkel-Poole, even when 
it is known that the rate of decay is proportional to more complicated 
voltage dependences than exp (V/V,). It is the purpose of this appendix 
to outline why this is so. For simplicity, we assume that conditions are 
such that only one decay mode is important. 

Usually we can write the decay of the stored charge, or, equivalently, 
that of the flatband voltage V, in the following form, 

dV 

Tt = ~U(V + bV b ) - f(0n (12) 

where, in turn, 

f{V) = Ae«< v '\ V = V + bV b , (13) 

A being a dimensional constant and g(V) a function of temperature 
and applied bias as well as of various powers of V (and possibly of 
log, V). The point of writing f(V) in the form given in (13) is to focus 
primary attention on the exponential nature of the functional depen- 
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dence of f(V) on V. Thus, f(V) can be approximated over a wider 
range of V by a Taylor expansion of g{V) in (13) than by a Taylor 
expansion of f(V) directly. 

Suppose at h « r re i B *, V = Vi. Then, assuming /(0) « f{V) for 
V of interest, we can readily integrate (12) if we approximate g(V) 
by g(V& + y(V — Vi), where 7 may be chosen in various ways 
(V't = Vi + bV b ). For example, we may desire upper and lower 
bounds for V(t), in which case we can use upper and lower bounds for 
g(V). In some cases, 7 = g'iYi) or 7 = \jg(V[) - ff(7a)]/(7i - V 2 ) 
give such bounds. In all cases, we obtain the characteristic decay of 
which eq. (3) is one example. 

A somewhat more appealing approach that avoids approximating 
quantities in exponents is the following. 28 ' 29 Again, assume /(0) « f(V) 
to write the integral of (12) in the form 



h h " lv x f(V) A J Vl 



vdVe-'Wg'iV) 



(14) 



g'(V) 

According to the intermediate-value theorem, this becomes 

-irm**"-*"' 91 (15) 

where V 2 < V 3 < Vi. (Physically, we must have g'(V) > 0.) Thus, 
given Vi and h, for each 7 2 of interest, U follows from (15) to within 
the uncertainty of choosing g'(Vs). (This will normally be, at most, a 
factor of 2 for V 2 ^ l/2Fi.) The envelope function V e (t) follows at 
once from (15) by setting h and exp [ m — g(Vl)'] to zero. 

t = e~«W/Ag'(yi). (16) 

If we note that V' e (t) = V e (t) + bV b , differentiating (16) with respect 
to V b at fixed t and assuming that g'{Vz) is only weakly dependent on 
Vb, it follows at once that 

dV e . 

Under these conditions, the vertical separation between envelope 
functions V e (t) corresponding to two different applied biases V\ and 
VI (at the same temperature) will be b{V\ — VI), even when these 
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envelope functions are not straight lines, but are given in general 
by (16). 

Equation (15) may be rewritten in the characteristic form as 

g(V 2 ) - g(Vi) = -log. [1 + Ae°Wg'{V*){U - <i)], (17) 

where we may write 

gW) - g(Vl)ttg'(Vi)(V 2 - 7i). (18) 

Equation (18) follows from the mean-value theorem (Va < V* < Vi). 
Thus, while it is clear that the exact solution will in general not be 
precisely given by (17) and (18), such a solution will still show the 
same general features as the above for physically reasonable g(V). It 
should be noted that, since g'{V' z ) and g'(Vl) actually represent aver- 
aged quantities, and since g(V[) is simply an initial condition, the 
results given in (17) and (18) are rather insensitive to the detailed 
dependence of g(V') on V. 

One final point should be noted. For V sufficiently small and Vb = 0, 
we may write 

^«-/'(0)7. (19) 

The solution to (19) is the common exponential decay {exp ( — t/t ), 
t = E/'(0)] -1 } now with voltage V an exponential rather than a 
logarithmic function of time t. For £» T rc i a x, clearly V(t) will be of 
this form. But this is well beyond the region in time of the primary 
decay of the charge, and it is this primary decay governed approxi- 
mately by (17) to (18), which is of crucial importance for evaluating 
charge storage devices. 

APPENDIX B 

B.I Decay through the Insulator 

Section II indicated that the electric field S in the insulator drives the 
charge-decay current [see eq. (1)]. This electric field arises from the 
stored charge and the applied bias. The purpose of this appendix is to 
express S in terms of two measurable voltages, the flatband voltage, 
which measures the stored charge, and the applied-bias voltage. 

Let Co be the capacitance of the oxide in series with the semicon- 
ductor, and let C,- be the capacitance of the insulator. We assume C 
(as well as C») is independent of voltage. (We thus ignore the bias 
dependence of the semiconductor capacitance.) Simple capacitative 
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division implies that, if V b is applied to the gate and there is no stored 
charge, then the voltage V b dropped across C, is given by 

V\ = V h yQ. (20) 

Vb Vh l/d + l/C K } 

On the other hand, if a charge Q is stored and V b = 0, then the voltage 
V a across C\ is given by 

V a = Q(Ct + C )-K (21) 

See Fig. lb. But Q is related to the flatband voltage V fb by the relation 

Q = CiVft. (22) 

See Fig. Id. Thus, 

V = Va, y^- (23) 

Va V/b l/Ci+l/C, K J 

The total voltage drop V T across C„ which produces S, will in general 
be a superposition of the two ; thus, 

V r -r.-V> 1/c %\ /Cm (V» + l>V*, (24) 

where 

b = Co/Ct. (25) 

In passing from eq. (1) to eq. (2) in the text, we have let V,(T) absorb 
the pref actor (24). Thus, V,(T) depends both on the physical processes 
of the decay and on the geometry of the device. 

It is important to note that b can be greater than as well as less than 
unity. In principle, b can be computed from knowledge of C„ and C t . 
It is more easily obtained from the measured decay curves, a more 
reliable method. Its slight temperature variation is somewhat of a 
puzzle. If at higher temperatures additional capacitances C\ and C 
arose in series with C, and C , respectively, then b would become b' 
given by 

One expects the semiconductor portion of the device C to be more 
susceptible to degradation than the insulator side C,-. If C' t = °c and 
Co < °°, then b' < b, as observed. For lower temperatures, b is close 
to its value predicted by (25). Fortunately, b is determined within 
experimental error at a given temperature to be independent of V b . 
Thus, at each temperature we can predict zero-bias behavior. We can 
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then predict room temperature, zero-bias behavior from higher tem- 
perature, and zero-bias behavior without knowledge of b. This is be- 
cause only the factor 67& enters eq. (2), and V& = for zero bias. 

8.2 Decay through the oxide 

If the primary discharging current is through the oxide (rather than 
through the insulator, as assumed throughout this paper), then the 
electric field that produces this current must be calculated from the 
voltage across the oxide. Using the notation introduced above, we 
note that the portion 7& of the bias voltage 7& dropped across the 
oxide is 

K-^' i/c'+i/c - (27) 

On the other hand, it is clear from Fig. lb that the contribution V' a of 
the voltage resulting from the stored charge is just Va given by (23). 
Thus, the total voltage drop V' T across C is given by 

r ' " f: ~ n - i/a + i/c. (V - n) - (28) 

Thus, unlike the case of Vt (24), no geometrical factor analogous to 
b enters. 

APPENDIX C 

In this appendix, we derive two minor results used in the text. 

(i) To show that Vjb (0 < V e (t), we need only compare the argument 
of the log e in eq. (3) with 1 1 = and V\ = 7J with that of eq. (4) when 
written 



7.(0 - V.log.[Yo + £eW-) H 



(29) 



From (3), we have 

7/6(0 = V 8 log e Mr-™* + i ^^^V l l ■ (30) 

As the argument of the log e in (30) is less than that in (29), it follows 
at once that 

7,»(0 < V e (t). (31) 

[Additional decay mechanisms important for very short times (<1 
minute), which are not included in (30), will, of course, reduce 7/6 
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even further.] It may be noted that Vi and t\ may be chosen so that 
V/b(t) given in (3) exceeds V e {t) and V fh approaches V e from above. 
Such unphysical behavior would result if one mistakenly chose a value 
of h too large ; that is, if, after completing the charging of the device 
(t = 0), one recorded Vi at h using a clock that read a time sufficiently 
larger than at t = 0. 

(it) To show that Vi « V\ when h « r re i a x, we calculate Vi from 
eq. (3), in which we set Vi = V\, h = 0, and V b = 0. Then it follows 
that 

*-'."*[•**•(» + £)"*]■ < 32 > 

If now ti « r re i ax (the usual case), then 

Vi w 7? - — F.. (33) 

7" re lax 

Since 7„ < Vi and h £^ 10 _2 r re i a x at most (for room temperature and 
zero bias 10~ 7 is typical), we may use Vi « V\ to a very good degree 
of approximation. The actual, initial, flatband voltage may be larger 
than Vi, owing to other decay mechanisms which may be important 
for t < 1 minute. It is Vi [or Vi obtained from Vi using eq. (3)], how- 
ever, with which we are concerned. 

APPENDIX D 
The "b-Factor" 

The disagreement between the predicted 6-f actor (25) of 2.0 and 
the much smaller measured 6-f actor (Table I) of 0.8 to 1.2, which we 
shall call b', calls into question our assumption that the electric field 
in the insulator is spatially constant. This follows because, if the field 
is spatially constant (and if the current per carrier depends only on 
the electric field, as is physically reasonable), then the measured and 
predicted 6-factors, b' and 6, would agree. Since, in fact, they do not 
agree, we conclude that the field is not spatially constant, a result, 
perhaps, of charge stored in the insulator. It is the purpose of this 
appendix to investigate some consequences of charge stored in the 
insulator (in addition to charge stored at the oxide-insulator interface) 
on the decay of the flatband voltage Vfb. 

Let us begin by asking what information our flatband-voltage mea- 
surements tell us in light of the possibility of having charge storage in 
the insulator. Referring to our empirical result (2), we note that, at a 
fixed temperature, the rate of decay of V fb (t) is a function of 
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CFf§(0 + b'Vb]. Thus, whatever electric field governs the decay of 
V/b(t), it too must depend on V /b and V b in the form \_V fb {t) + &'F&]. 
The next step is to express the various fields and voltages of interest 
in terms of the stored charge. Let Q,(t) be the charge stored at the 
oxide-insulator interface at time t, p f (x, t) the density of charge stored 
in the insulator, e, the dielectric constant of the insulator, and e that 
of the oxide. It follows (see Fig. Id) that the flatband voltage is 
given by 

V /b (t) = W- + r ^^ fa - x)dx (34) 

W JO *i 

- V,{t) + V m {t) (35) 

' di PT(x,t) 



i; 



(di - x)dx (36) 

= V T (t), (37) 



where 



Pt = P.+ (Q,/A)S(x). (38) 

Here A is the cross-sectional area of the device, d, is the thickness of 
the insulator, 8 (x) is the usual delta function, V, is the contribution of 
Q a to Vf b , and V m (the w referring to moment of charge) is the con- 
tribution of the charge stored in the insulator to V fb . (Throughout our 
discussion, we are assuming that no charge is stored in the oxide.) A 
straightforward but more involved calculation leads to the electric 
field in the insulator, < x < di, under general charging and applied- 
bias conditions : 

E(x, t) - (1 + b)-[ ?iWM + b JL» + j- dx , fJ ££ ( x + 6 1) 

If we integrate (39) or (40) over the insulator, we obtain the average 
field times d„ the voltage drop across the insulator V d {t) given by 

V d (t) = (1 + 6)->[F„(0 + F m (i) + bV b -] 

= (1 + 6)-'[F /6 (0 + 67 J. (41) 
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Using these results, we can discuss what physical effects can and cannot 
lead to the observed behavior (b' < b). 

Owing to the presence of charge in the insulator, the electric field 
at the oxide-insulator interface will be reduced from its value in the 
absence of such charge. If 2?(0+, t) controls the decay, perhaps this 
reduction can lead to a reduced b-factor. Unfortunately, this effect 
would result in (&' > b) rather than {b' <b). To see this, we note 
from (39) that 

E(0+, t) = (1 + 6)- 1 df 1 IT.(0 + bV b - bV m (t)~] 

= T& i [ b b v ^- b ' VM+b ' Vb } (42) 

If this expression is to be a function of [7/6 (0 + b'V h ~], it is necessary 
that 

V fb (t) = (6'/6)7.(0 - b'V m (t). (43) 

To satisfy both (35) and (43), it is necessary to have 

"-(0 = 1/ 1 6 / ~ + 1 f vmo m 

and 

r 'W-17RTT 7 ' iW ' (45) 

However, if b' < b, then 1/6 < 1/6', and the coefficient of V /b (t) in 
(44) is negative. This is rather disturbing since, as a result of the initial 
charging and subsequent decay of the charge through the insulator, 
one would expect V m (t), as defined in (35), to be positive. In addition, 
(44) implies that, as V/i decreases during charge decay, V m must 
increase algebraically. This can occur only if charge stored at the 
oxide-insulator interface becomes trapped in the insulator neutralizing 
the charge trapped there. This effect would be enhanced if more charge 
were neutralized near the interface (x = 0) than in the middle of the 
insulator. Although this may be the source of (6' < 6), the origin seems 
physically unreasonable of the bulk trapped charge of sign opposite 
the stored charge which, although not neutralized by the charging 
current, is neutralized by the decay current. We conclude that, if 
(6' < 6), then the field at the interface probably does not control the 
charge decay. 

We note in passing that, since the voltage drop V d {t) across the 
insulator (41) is already a function of [V fb {t) + 67&], it also cannot 
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be the average field in the insulator which is controlling the decay. 
This hardly requires further elaboration. 

One effect that charge stored in the insulator has on the decay is to 
enhance the electric field near gate (x = di) relative to the field near 
the interface (x = 0). This enhancement may be sufficient to remove 
charge stored in the region x < x < di of the insulator during the 
initial portion of the decay, leaving significant stored charge only in 
the region < x < x of the insulator and at the interface. If, during 
the remainder of the decay, the decay rate is governed by the electric 
field at x , we then expect to observe (&' < b). Let us see how this 
comes about. 

If pi(x, t) = for x > x , then the electric field at x = x is, accord- 
ing to (39), given by 



where 



+ 6 _[( 6 + l)^_ 67m(0 ]|, (46 ) 
Q t (t) = A f di dx'piix', t) (47) 

Jo 

= A f XO dx' Pi (x', t), (47a) 

Jo 

the size of the stored charge in the insulator. Alternatively, using (40) 
for the field, we obtain 

*<*« () = (TT5J3; [ - VVt{1) + ^ (1 + b) ™ + H ■ (48) 

where Q(t) is the total stored charge defined by 

Q(t) =Q.(0 + Q,(0. (49) 

For E(x , t) to be a function of (V fb + b'V b ), it is necessary that 
[using (37)] 

V T = -b'V T + £ (1 + b) Q& (50) 



or that 



v (t) - 1 + 1/b Q(t) m\ 

VT{t) -\TTJb'~c7' (51) 
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Now (&' < 6) implies that C,-7r(<) < Q(t), which means in turn that 
some charge is stored in the insulator, a consistent result. 

Having seen how a measured 6-factor (6') can arise which is less than 
the computed 6-factor (6), we must inquire as to whether the effect is 
sufficient to explain the measured results. For simplicity, let us assume 
that pi is uniform for < x < x„. Then it follows that 

Vt = % + % (1 " Xo/2di) - (52) 



(53) 



It then follows from (48), (51), and (52) that 

Q/d _ 1 + Qt/Q. = 1 + W 

V T 1 + «2t/Q.)(l -x./2dt) 1 + 1/6 

If we consider the case at hand where 6 = 2.0 and 6' is essentially 1, 
then (53) implies that x must be at least 0.50 X d i} that is, that the 
stored charge in the insulator must extend at least 50 percent of the 
distance from the interface to the gate. If Qt/Q, is 1.0, then x = di. 
Thus, Qt/Q, must be at least 1.0 for, if it were smaller, then x„ > d if 
which is not possible. Therefore, if the charge stored in the insulator 
is uniform between x = and x = x , then it is necessary to under- 
stand the observed 6' that 0.50 < x„/di < 1 and 1.0 < Qt/Q», the 
latter implying that more charge must be stored in the bulk at the 
insulator than at the interface. 

In the preceding paragraph, we have assumed that the stored charge 
was uniformly distributed in the region < x < x . In fact, we expect 
the charge to be more dense near the interface (x = 0) where the field 
is lowest than near x = x . To satisfy (51), this would require larger 
values of Qt/Q, and of x than for the uniform case. We have also 
assumed that it is the electric field at x that controls the decay. It is 
possible that E(x, t) for x < x in fact performs this function. This 
would further increase the values of Qt/Q, and x required to achieve 
(6' < 6). One's lattitude here is rather limited, however, for, as we 
have seen, if it is E(0 + , t) which controls the decay, then Q t /Q, becomes 
negative. We offer the above, therefore, as possibilities only. 

Another source of the (6' < 6) effect is that the charge may be 
extraction-limited, that is, controlled by the field at x = d { . If we put 
x = di in eqs. (46) and (47a), then we again obtain (51) relating the 
total stored charge to its first moment. We noted above that we could 
understand the measured 6-factor for x = di and a uniform, stored, 
insulator charge if Q t /Q, = 1.0, a reasonable but perhaps somewhat 
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large value. However, if the current is extraction-limited, then we 
expect that the charge would be more dense near x = di than in the 
bulk of the insulator. As a hypothetical example, suppose that the 
stored, insulator charge is uniform for x\ < x < di and zero for < x 
< Xi. Then from (37) we obtain 

v _Q. , Q.1A x x \ .... 



(55) 



It then follows from (49) and (51) that 

Q/Cj = 1 + Q t /Qs 1 + 1/6' 

V T " 1 + (Q,/Q.)(l - *i/d«)/2 1 + 1/6 ' 

This relation provides a much greater possibility for obtaining (6' < 6) 
than (53). For example, if Q t » Q s , then 

6' S (1 + 2/6)- 1 < 6/2. (56) 

Or, if Zi P^ d„ then 

6' - 6[(1 + b)Qt/Q. + l]- 1 , (57) 

from which 6' = 1 would follow for 6 = 2.0 if Q,/Q B = 0.333, a very 
reasonable value. For xi < c/„ somewhat larger Q t /Q s would be re- 
quired to satisfy (55). However, in most cases Qt/Qs would be smaller 
than that in the previous example (53), in which the stored charge 
was assumed near the oxide-insulator interface. We conclude that, 
while we have indicated the possibility of how (6' < 6) can come about, 
further work is required to really pin down and calculate the measured 
6-factor 6'. 
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